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Abstract: CP violation in the standard model originates from the Cabibbo-Kobayashi- 
Maskawa mixing matrix. Upon integrating all fermions out of the theory, its effects are 
captured by a series of effective nonrenormalizable operators for the bosonic gauge and 
Higgs fields. We compute the CP-violating part of the effective action to the leading 
nontrivial, sixth order in the covariant gradient expansion as a function of temperature. In 
the limit of zero temperature, our result addresses the discrepancy between two independent 
calculations existing in the literature [1, 2]. We find that CP violation in the standard model 
is strongly suppressed at high temperature, but that at T < 1 GeV it may be relevant for 
certain scenarios of baryogenesis. We also identify a selected class of operators at the next, 
eighth order and discuss the convergence of the covariant gradient expansion. 
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1 Introduction 

The combination of the discrete parity and charge conjugation transformations (CP) is 
a symmetry of the gauge interactions of the standard model. However, a complex phase 
of the mixing matrix appearing in the fermion-scalar interactions breaks this symmetry. 
Combined with the breaking of baryon number conservation through the chiral anomaly 
[3], this opens up the possibility of explaining the observed asymmetry between matter 
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and antimatter in the universe, using only known electroweak-scale physics. This scenario 
is commonly referred to as electroweak baryogenesis [4]. 

The dynamical processes leading to permanent change of baryon and lepton numbers 
in the standard model are intrinsically nonperturbative. As a consequence, we cannot ex- 
pect to compute electroweak baryon number creation rates within a purely perturbative 
scheme. One approach to deal with this problem is to evolve the entire system numerically, 
for instance on a spacetime lattice. The advantage of this is that the out-of-equilibrium 
dynamics would be correctly treated, but a serious complication arises from the fact that 
the fermions need to obey quantum dynamics. This is numerically very challenging, as in 
addition to the usual issues associated with lattice fermions one needs to evolve each quan- 
tum fermion mode separately in time [5-7]. Fortunately, the gauge dynamics responsible 
for electroweak baryon number violation is very well understood numerically, both in and 
out of equilibrium, and can be convincingly described by classical dynamics or stochastic 
equations. The reason for this is that a change in baryon number is an infrared process, 
mediated by "spatially large" gauge field configurations like sphalerons. 

For the reasons listed above, it is natural to regard CP violation as a small, pertur- 
bative bias originating from the back-reaction of quantum fermions, living in a classical, 
nonperturbative gauge field background. Consequently, one may integrate out the quantum 
fermions in the path-integral formulation of the standard model, creating a purely bosonic 
effective theory of gauge and scalar fields. The effects of CP violation are then captured 
in a series of effective, nonrenormalizable operators. These can in turn be computed using 
a perturbative (or otherwise appropriate) expansion, their form and coefficients becoming 
functions of the fermion masses and mixing angles. 

The benchmark calculation of the effective CP violation in the standard model was 
performed by Shaposhnikov [8]. Rather than integrating out fermions, his aim was to 
calculate the coefficient Scp of the leading CP-violating operator, defined in eq. (4.1), in 
perturbation theory in terms of diagrams with both bosonic and fermionic internal lines. 
The conclusion was that at asymptotically high temperatures 6cp — JA/T^^, while at zero 
temperature Jcp — JA/v^"^, where v ~ 246 GeV is the Higgs field vacuum expectation 
value and A denotes a specific combination of fermion masses, 



for quarks and analogously for leptons. In addition, J stands here for the Jarlskog invariant, 
defined in terms of the 3x3 fermion mixing matrix V (repeated indices not to be summed) 



where Cij = ^k^ijk- In the quark sector, V represents the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix, and J ~ 2.96 x 10~^ [9]. Inserting the quark masses according to ref. [9], 
rriu = 2.3 MeV, nid = 4.8 MeV, = 1.275 GeV, m, = 95 MeV, nit = 173.5 GeV, 
mfe = 4.18 GeV, one obtains (5cp ^ 10"^*^ and lO'^^ for T = 100 GeV and T = 0, 
respectively. In combination with later work [10-13], this simple estimate is the origin of 
the common lore that standard model CP violation is not sufficient to explain the observed 



(1.1) 





- 2 - 



baryon asymmetry, usually phrased as the net baryon-to-photon number ratio 

^^6x10-1°. (1.3) 

Some time ago, Smit [14] demonstrated that the situation is in fact somewhat more 
complicated in that at least at zero temperature, integrating out the fermions does not 
lead to CP-violating operators suppressed by powers of the small Yukawa couplings, but 
that instead 5cp — J ■ The first attempts to explicitly evaluate the leading CP-violating 
operators and their couplings, however, resulted in controversy: the operators found by 
Hernandez et al. [1] break P and conserve C, while all operators reported by Garcia- Recio 
and Salcedo [2] conserve P and break C. 

While the non-suppression of effective CP violation at zero temperature is highly en- 
couraging for baryogenesis, it is important to recall that realistic scenarios do not take 
place at zero temperature. "Hot" electroweak baryogenesis typically operates at tempera- 
tures around 100 GeV [4, 15], while "cold" electroweak baryogenesis in an infrared-heavy 
out-of-equilibrium environment is assumed to take place at temperatures below roughly 
40 GeV [16-19]. (The temperature prior to electroweak symmetry breaking may be ar- 
bitrarily low, if emerging directly from an inflationary stage.) Hence, knowledge of the 
temperature dependence of effective CP violation in the bosonic sector of the standard 
model would clearly be desirable. To this end, our present work is aimed at generalizing 
the calculations of refs. [1, 2, 14] to nonzero temperature and discussing their implica- 
tions for the cold electroweak baryogenesis scenario. Some of the results derived here were 
presented already in ref. [20]. 

The outline of our paper is as follows. In order not to obscure the physical picture with 
unnecessary technical details, we begin by providing a qualitative sketch of the calculation 
in section 2. Section 3 on the other hand presents our results in detail; in particular, 
section 3.1 reviews the temperature dependence of the leading CP-violating operators as 
reported in our previous paper [20], while section 3.2 provides a partial list of subleading 
operators at zero temperature, and section 3.3 investigates the effect of an additional source 
of CP violation in the lepton sector. The implications of our findings for baryogenesis are 
finally discussed in section 4, while the technical details of our computation are relegated 
to several appendices. For the benefit of the reader, we also briefly summarize our notation 
and conventions in appendix A. 

2 Flow of the argument 

As explained above, our aim is to integrate out the fermions from the full standard model 
in order to derive an effective action for the gauge and Higgs fields, and finally to identify 
the CP-violating part thereof. As the effective action cannot be determined in a closed 
form for arbitrary spacetime-dependent background fields, we will resort to a covariant 
gradient expansion in its evaluation. As was argued already in refs. [1, 2], this is expected 
to provide a good approximation for temperatures not exceeding the charm quark mass. 
The results presented in section 3.2 will in addition allow us to directly investigate the 
convergence of the expansion. 
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2.1 Gradient expansion of the effective action 

Let us start by considering a generic Euclidean field theory of chiral fermions, denoted 
collectively as ipix) = {ipi^{x) , ijjji{x)) , coupled to external gauge as well as scalar fields. 
The partition function of such a theory can be expressed as a functional integral over the 
fermionic fields, 



Z 



[d^] [dtp] exp 



ip{x)D{x)ip{x) 



(2.1) 



where the effective action F is formally given in terms of the Dirac operator D{x) as 
r = — Trlog-D. However, as shown in refs. [21, 22], its evaluation is subtle due to renor- 
malization ambiguities and the contribution from the chiral anomaly. To this end, we write 
the Dirac operator in the chiral (left-right) basis in a block form. 



D 



-0L "ILR 



(2.2) 



and define the object K = mLRfn-RL — -0L"iRL-^R"^RL- The effective action is then given 
by the sum of a normal and an abnormal parity component, and F", expressed in terms 
of -ff as 

r+ = -lReTr(log/C), T" = Im 1^(75 log if) + Tgwzw, (2.3) 

where Tgwzw is the anomalous Wess-Zumino-Witten term [21, 22]. 

The dominant source of CP violation in the standard model is the CKM matrix in the 
quark sector, to which we will restrict our attention from now on. (The effect of additional 
CP violation originating from the lepton sector will be discussed in section 3.3.) In the 
block form indicated in eq. (2.2), the quark Dirac operator reads 

(iI)u + $ + Mi^ 



D 



\ 



V " 

















\ 



IPu + 
lPd + 0j 



(2.4) 



cf. eq. (A. 2). Here, Z^, VF^ and are the weak intermediate boson and gluon fields, 
respectively, (p is the Higgs field, and M^^d stand for the complex mass matrices of u-type 
and d-type quarks before diagonalization. Further notation, including explicit expressions 



for the covariant derivatives D^'*^, is summarized in appendix A. Throughout this paper, 
we use the unitary gauge, in which only a single physical scalar field (j) appears. 

Using the explicit expression for the quark operator (2.4), we obtain for the operator K 



K 



\4>''/v^)MuMl-{Ip^ + $){Ip^ + ^) -T^+(0^ + yj) 

-ir"(0„ + ^) {<t>''/v'')MdMl-{Ipd-$){U>d + i>), 



(2.5) 



with ipfj_ = (p ^[df^, (j)]} The gluon fields have now been omitted for the sake of simplicity, 
as they only contribute at the subleading, eighth order in the gradient expansion [23] . One 



^Recall that at this stage of the computation, we are deahng witli (pseudo-)difrerential operators. All 
fields act as multiplicative operators and the commutator [d^ , <j}] thus corresponds to a multiplication by a 
derivative of the Higgs field. 
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nevertheless has to bear in mind that the trace in eq. (2.3) still involves the color space, 
thus giving rise to the trivial prefactor Nc, the number of colors. 

In the next step one expands the logarithm of K in powers of its off-diagonal part, or 
equivalently, in powers of the fields. Following ref. [2], let us introduce a shorthand 
notation for the inverse of the diagonal components of K, 

= ^^M^id - i^u4 ± ^)i^u4 + N-} = ^M^^Kd+P'- (2-6) 

One obtains a compact expression {logK) = '^'^^oi^og K)2m where (O) stands for TrO 
or Tr(75C') in F"*" and T~ , respectively, and 

(log A02n = ([iV^Tf +(0rf + ^)7VdTr"(0„ + y^)]") . (2.7) 

The trace is evaluated with the help of the method of covariant symbols, explained in 
detail in appendix B. In the covariant gradient expansion, the result is further expanded 
in powers of the covariant derivatives and the Z^, ip^ fields. All the contributions are then 
classified as (log K)2n+ni where, as above, 2n is the number of fields and m counts 
the covariant derivatives, and ip^i fields. The order of a given operator in the gradient 
expansion is 2n + m. 

In order to have any CP violation, at least four fields are needed. This is because 
the simplest CP- violating observable built from the CKM matrix is the Jarlskog invariant, 
eq. (1.2), which is composed of a product of four CKM matrix elements. Therefore, CP 
violation can only appear at order four or higher in the gradient expansion. However, 
as was shown by Smit [14], order four actually gives no CP violation, a result which we 
verified to hold also at nonzero temperature. Since the chiral anomaly only contributes at 
this order in the gradient expansion, it can be discarded as long as only the CP-violating 
part of the effective action is desired [2]. We can thus use eq. (2.7) to identify all effective 
bosonic CP-violating operators of the standard model. 

At zero temperature (and in an even-dimensional spacetime), only operators of even 
orders in gradients appear in the action as a consequence of Lorentz invariance. Nonzero 
temperature can in principle lead to operators of the 4+1 type, but we checked by an 
explicit computation that they do not contribute to CP violation in the standard model. 
The leading order at which CP violation appears is then order six. In fact, there are no 
contributions of the 2n+0 type for any n, as demonstrated by Garcia-Recio and Salcedo [2] 
and verified by us also at nonzero temperature. Therefore, all contributions at the leading 
nontrivial order must be of the type 4+2; they are investigated in detail in section 3.1. In 
order to assess the convergence of the gradient expansion, we evaluated a selected class of 
operators at the next, eighth order at zero temperature. There is no CP violation coming 
from operators of the type 8+0, which leaves us with two possibilities: 6+2 and 4+4 type 
operators. In section 3.2 and the associated appendix D, we derive the full CP-violating 
action in the 6+2 sector. We leave aside the evaluation of the 4+4 operators since the 
immense number of possible combinations of Lorentz indices and fields makes it impossible 
to even list all contributing operators in a paper. We merely note that the first CP-violating 
and P-violating operators calculated in ref. [23] are of the 4+4 type. 
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2.2 Remarks on the gradient expansion at nonzero temperature 

While the above discussion lays down a concrete path both at zero and nonzero tempera- 
ture, there is one additional subtlety in the calculation that needs to be addressed. It is 
namely a well-known fact that the gradient expansion breaks down at nonzero tempera- 
ture. This may be naively understood as a consequence of the periodic boundary condition 
imposed on the quantum fields in the Matsubara formalism; for any given bosonic field 
(j) as a function of the imaginary time r, this reads (j){T = 1/T) = (j){T = 0), where T 
is the temperature. Expanding (^(r = 1/T) — (j){T = 0) in Taylor series around r = 0, 
we obtain a series of terms of successively higher order in the gradients, which only pro- 
duce zero upon the full resummation. Truncating the expansion at any finite order gives 
spurious contributions that, apart from the periodicity itself, may break other symmetries 
of the effective action. The fully gauge invariant effective action is necessarily nonlocal; 
the compactification of the time dimension results in the appearance of the Polyakov loop 
operators [24, 25]. 

When computing the bosonic effective action of the standard model at nonzero tem- 
perature, the above problem shows up in those terms containing temporal derivatives of 
the various fields. While the proper way to deal with them would be to follow the strategy 
of ref. [25], we took a simpler approach, avoiding the necessity to cope with the unphysical 
contributions that possibly break the symmetries of the effective action. Concretely, we 
performed on purpose the gradient expansion also in the forbidden temporal derivatives 
and restricted our attention to the Lorentz-invariant part of the result, which does not suf- 
fer from the above problem. It can be defined unambiguously by discarding all operators 
containing temporal indices (whether in derivatives or gauge fields), and then replacing the 
remaining spatial indices with Lorentz spacetime ones. This should be kept in mind when 
inspecting the results of section 3.1. 

2.3 Momentum integration 

The evaluation of the effective action using the method of (covariant) symbols leads to 
momentum integrals of the type 



where the integration symbol represents either an integral over all spacetime components 
of momentum (at zero temperature) or an integral over spatial components together with a 
Matsubara sum over the temporal component (at nonzero temperature). For full details of 
our notation, see appendix A. The dependence of the result on the Higgs field (p is trivial, 
as can be seen by rescaling the spatial momentum phy (j)/v and redefining the temperature 
to Tefj = Tv/(f). We thus see that the above momentum integral can be written in the form 
{(f) / v)^ where \ = d + 2{k + 1) — 2 J21^=ii''^i + is the spacetime dimension, 




p 
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and 



"^1 ) ■ ••I'^rtt^'n, 



Jp 



(2. 



X tr 



The temporal component of momentum in this integral is determined by Tefj rather than 
T itself; this should be kept in mind throughout the detailed derivation of some of the 
thermal integrals in appendix C. An analogous argument applies at zero temperature, 
the only difference being that the whole four-momentum p must be rescaled by (j)/v. In 
addition, Lorentz invariance then ensures that the integrals Irfti tn ^PPsar only in 
specific combinations, discussed at length in section 5 of ref. [2], 



1 



= i:C)4'i-;....„,.„. (2.10) 

In the effective action, this is accompanied by the factor {(f) /vy where now \ = d + 2k — 

As the CP transformation is equivalent to complex conjugation of the mass matrices, 
Mu,d, the CP-violating part of the effective action is proportional to the imaginary part of 
the momentum integrals, 7^/., . = ilm75/+, ^ + and analogously for l!^^ x ^ + . 
At zero temperature, the leading-order CP-violating part of the effective action (belonging 
to the 4+2 sector) is proportional to the single integral if 12 2- "^^^ action in the 6+2 
sector turns out to be proportional to a single integral as well, this time 7^i 1 1 2 21 some 
details of its calculation are presented in appendix D. The evaluation of the temperature 
dependence of the leading-order result on the other hand requires the determination of the 
integrals 7{^'i 2 2 ^^^^ -^1123 -^1132' with non-negative integer superscripts. Their 
calculation is detailed in appendix C. 



3 Results 

In this section, we list all of the results we have obtained at different orders of the gradi- 
ent expansion both at zero and nonzero temperature. In practice, the computations were 
performed using the Mathematica package Feyncalc [26] . Our code uses the method of co- 
variant symbols to evaluate the trace in eq. (2.7), taking as input the expansion coefficients 
from eqs. (B.8) and (B.15). The resulting momentum integrals were simplified manually 
and subsequently evaluated numerically. While the general derivations presented in this 
paper are formulated and valid in a Euclidean spacetime of an arbitrary dimension, the 
results presented in this section are specific to the physical case of four dimensions. 

It is worthwhile to remark that in order to cast the effective action into the canonical 
form presented below, integration by parts (in coordinate space) was sometimes necessary 
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to avoid, for instance, the appearance of operators with all derivatives acting on a single 
field. When doing so, we tacitly assumed that no surface terms arise from such integration. 
This is equivalent to demanding that the fields are topologically trivial, which is consistent 
with the fact that we are performing a gradient expansion. On a more technical note, when 
integrating by parts, one has to pay attention to the fact that the mass matrices in the 
propagators (2.6) contain a factor of the Higgs field (p. At zero temperature, this can be 
easily gotten rid of by rescaling the momentum variable prior to the integration. This is no 
longer possible at nonzero temperature, and we carried out the integration by parts before 
the step leading from eq. (2.8) to eq. (2.9). This results in particular in the appearance of 
momentum integrals of the type (C.5). 

3.1 Leading order 

At zero temperature, the CP-violating part of the effective action was evaluated to its 
leading, sixth order in refs. [1, 2]. The main goal of our preceding paper [20] was to 
resolve the discrepancy between these two calculations and to generalize the results to finite 
temperature. As observed in the previous section, all relevant CP-violating operators at 
this order are of the 4-1-2 type, and hence the resulting (Euclidean) effective action takes 
the form 



and the parameter A was introduced in eq. (1.1). 

The operators On appearing above are composed of the W^, and (p^ fields; the 
subscript n counts the number of the latter two. They can be further split into P-even and 
P-odd parts. On = O^ + 0~ . At nonzero temperature, the effective action is no longer 
Lorentz-invariant due to the presence of the thermal bath, but can be formally written in 
a covariant way, if one introduces a timelike vector defining the rest frame of the thermal 
bath, = 5fj,o. We can then also divide each operator On into its Lorentz-invariant part, 
containing no u^'s, and the rest. We have verified that the Lorentz-noninvariant part of 
the action vanishes in the limit of zero temperature, as it must. The Lorentz-invariant part 

^We write the integration measure explicitly to emphasize that Kqp^ is defined at zero temperature; the 
temperature dependence of the CP violation effects originates through the coefficients Ci introduced below. 
The coefficient is expressed in terms of a momentum integral of the type (2.10) as iJGpK^^^ = /i,i,2,2- 
Due to the hierarchy of the quark masses, it can be well approximated as K^p^ l/(167r^GFm^) ~ 334. 





where Gf = 1/(^2^^) is the Fermi coupling, the coefficient defined as 




(3.2) 
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Teff/GeV 

Figure 1. The coefficients C1-C13 plotted as functions of the effective temperature Tcff ~ Tv/(t). 
The thick lines correspond to the smallest and the largest of those q's that approach one at zero 
temperature, that is, ci and ciq. This figure was previously published in rcf. [20]. 



is on the other hand explicitly given by Oq = = C2 ~ ^ ^^'^ 

= - j{w+fw~^w~, + ^-^{w^fw~,w^, - 'j{w+fw-,w-^ 

+ ^W;W+W-Xa - ^W+W+W-^W~^ (3.3) 

- 2c,W;W^W~^W-, + ^W+W+W-Xu - c.c, 

Ot =1{Z^ + iff.) [c^{W+fW;W~, - c^{W+fW~W-, - c^{W+fW-W-^ 

- c^{W+ ■ W-)W;W-, + c,{W^ ■ W-)W^W-, + c,W+W^W-W-, (3.4) 

- ci2(W^+ • W-)W+W~^ - c^2W+W+W-W-^ + ci^W-W+W+W-,] - c.c, 

Ot =A{Z^Z, + ^^^,) [c^{W+fW-W- - c^{W-fw;W+] 

16 4 
--{Z-^) [c^{W+ -W-f- 2cQ{W+)\W-f] + -{Z^^, + Z,^^) (3.5) 

X [ci^{W+fW;W- + ci^{W-fW^W^ - 2cii(VF+ • W-){W+W- + W+W;)] , 

where "c.c." stands for complex conjugation. The coefficients Cj are complicated functions 
of the quark masses, temperature and the local Higgs field (i){x), of which the latter two 
appear in the combination Tcff = Tv/(j). The list of explicit expressions for the coefficients 
is lengthy and we thus provide it at the end of appendix C, in eq. (C.14). In the zero- 
temperature limit, the coefficients ci-cn approach unity while C12, C13 vanish, reducing our 
result to that of ref. [2]. We emphasize that all the operators written above are P-even and 
C-odd, generalizing the observation of ref. [2] to finite temperature. 

The behavior of the coefficients Cj as functions of Tefj is displayed in figs. 1 and 2, 
reproduced from ref. [20]. Our main observation concerns the rate at which these functions 
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10-2 r 
10-4 r 

10-6 r 

10-8 r 
10-14 

0.0 0.5 1.0 1.5 2.0 2.5 

Figure 2. The dependence of the c^'s on the Higgs field 0, plotted for three different values of the 
temperature, T = 1/3 GcV, 1 GcV and 3 GcV. In each case, the grey band is spanned by ci and 
cio- This figure was previously published in ref. [20]. 

fall as one progresses towards higher temperatures: while at Tcfj < 100 MeV the coefficients 
display a relatively slow variation, at phenomenologically interesting higher temperatures 
they start to decrease very rapidly with Tgff. The largest coefficient at all temperatures 
is cio, which reaches its maximal value of approximately 1.3 at T^s ~ 18 MeV, while at 
1 GeV it has already decreased down to the 10~^ level. At even higher temperatures, we 
observe a temperature dependence of the coefficients which is compatible with the estimate 
of ref. [8]. 

As highlighted in figure 2, the dependence of the coefficients q on the parameter 
Tcff has the impact that at any fixed nonzero temperature, their magnitude depends very 
strongly on the value of the Higgs field (j). Small values of the Higgs field correspond to 
stronger suppression of CP violation, while those regions of space where (j) is large witness 
an enhancement of CP violation. As discussed at some length in ref. [20], this observation 
has the effect of improving the convergence of the gradient expansion, as it alleviates the 
problems numerical simulations of standard model dynamics face in those regions of space, 
where the Higgs field changes its sign. 

3.2 Next-to-leading order 

The motivation to investigate the next-to-leading-order terms in the gradient expansion 
is twofold. First, as observed in ref. [2] and confirmed by us, all CP-violating operators 
at the leading order conserve parity. This is a very nontrivial and unexpected result, 
which is at odds with the findings of ref. [1]. At present, we cannot offer a satisfactory 
explanation for why all P-odd contributions, though present in the intermediate stages of 
the calculation, should cancel in the final result. It is then important to try to detect the 
first P-odd operators in the gradient expansion. A selected class of such operators, of the 
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type (W^W~)'^Z'^D, was indeed found by Salcedo at the eighth order [23]. His result, 
albeit incomplete, is certainly valuable, calling for both verification and generalization. A 
second, altogether separate, motivation for the next-to-leading-order calculation is that it 
allows us to examine the convergence of the gradient expansion. This is clearly a highly 
relevant issue for all practical applications of our results. 

As observed in the previous section, the order-eight operators are either of the type 
6+2 or 4+4, the former arising from the n = 3 term of eq. (2.7). At zero temperature, the 
full result in the 6+2 sector reads 



4 



r^t-^odd = -^NJl,^i,i,2,2 [-^ ) (Oo + Oi + O2) , (3.6) 

where again Oq = = O2 = and 

0+ =^W+f{W+ ■ W-){W-^W- - AW-W- + W-,W-^) + 2{W+fW+W- 

X (W~ W~ + W~ W~ + W~ W~ - AW~ W~ + W~ W~ + W~ W~ ) f3 7) 
+ 6(VF+ • W-)W+W+{-W-,W-^ + '^W-^W-^ + 2W-^W-, - W'Xu) 
+ AW+W+W+W^W-^W-f, - 6W;W+W+W^W-,W^^ - C.C., 

ot =2{z^ + [-2iw+fiw-fw;w- + 3{w+f{w-)^w;^w- 

- 2{W+)\W-fW+W-^ - A{W+)\W+ ■ W-)W-W~, 

+ Q{w+f{w+ ■ w-)w-w-, + (H^+)2(Ty+ . w-)w-w-^ 



+ 6(W^+ • W-fw;W-, - 9(Ty+ • W-fWtW-, + (t^+ • W-fW^W-^ (3.8) 

+ {w^fw;w-w-w-^ + {w^fw-w^w-w-^ + {w^?w~w-wtw-^ 

+ i^W-fW^W^W^W-^ - 3(T^+ • W-)W;w:jW-W-^ 

- 3(W^+ • W-)W;W-WtW-^ + 2(Ty+ • W-)W-WtWtW-^ - C.C., 

ot =IQ{Z^Z, + ^^^.){w+ ■ W-) [{w-fw+w+ - {w+fw~w-] 

+ 32(Z • (^)(VF+ • W~) [{W^ ■ W-f - {W-^f{W~f] (3.9) 

+ 2{z^ip, + z,^^){-{w+f{w+ ■ w-)w-w- - {w-fiw+ • w-)w;w+ 

+ [A{W+ ■ W-f - 3{W+f{W-f] {W+W- + W+Wp}. 

Note that similarly to the type 4+2 operators, there are no P-odd contributions. The 
prefactor of the result is given by a new type of integral, the detailed derivation of which 
is performed in appendix D. Similarly to the prefactor in eq. (3.1), it can be expressed in 
terms of a dimensionless parameter, K^p^, defined as (4/15)lf 2,2 = iJGpK^^. We find 



^6+2 



4 A r d^p ^ ^^2 65 . 10^ (310) 

15 G| J (2vr)4 n^^,(p2 + ^2)2 2^ (^2 + ^2_^)(^2 + ^1^.) ' 



see eq. (D.5). Interestingly, this integral is strongly infrared sensitive. It would diverge if 
we naively tried to set the u- and d-quark masses to zero, and in fact, 99% of its value comes 
from the i = j = 1 term in the sum. As in the case of Kqj?, a very good approximation is 
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obtained by factorizing out all the heavy quark masses, producing 

6+2 iKdP A rus 197\ g 

Kpp ;r— 75 7: log W 2.64 X 10 , (3.11) 

307r2G|m2m2 V m 120y ' ^ ^ 

where m = (m„ + md)/2. 

We have not determined the exhaustive list of CP-violating operators of the 4+4 type. 
The reason for not having completed this calculation is twofold. First, the result is plagued 
with ambiguities due to the existence of relations among different operators arising from 
integration by parts. (The approach used in ref. [23] seems more suitable for this purpose as 
it derives the effective action from the conserved current and thus avoids such ambiguities 
by construction.) This is, however, merely a technical obstacle, which can be overcome with 
some effort. The reason we have not attempted to do so is that the full CP-violating action 
in the 4+4 sector contains several thousand independent operators. While the coefficients 
of some of them probably vanish due to accidental cancellations, the result would still be 
of little practical use for subsequent numerical simulations, and would almost certainly not 
be possible to publish in a journal article. 

Finally, let us remark that at order eight, the gluon fields start to contribute to the CP- 
violating operators. All terms including gluons have been identified by Salcedo [23], and we 
have confirmed his result. The possible form of such operators is strongly constrained by 
symmetries: two factors of the gluon field strength tensor are needed in order to construct 
a color singlet, which together with the four fields necessary for CP violation make 
for an order-eight operator; no or (p^ fields can thus appear. The final result for the 
gluonic part of the Euclidean CP-violating effective action reads 

Tcp'oZ = ^A'i,2,2 (^)' [{W^?W;W- - {W-fw;:W;^] tr(G„,G„.), (3.12) 

where G^i^ = [(9^ + G^, dy + Gy] is the gluon field strength tensor and the trace is performed 
only over the color space. Our result is in agreement with that of ref. [23]. 

In order to investigate the convergence of the gradient expansion both at zero and 
nonzero temperature, we use the result for the effective coupling in the 6+2 sector (3.10), 
and compare it naively to K^p^, eq. (3.2). Generalizing both expressions to finite temper- 
ature by simply replacing the four-dimensional momentum integrals by sum-integrals, we 
define a parameter with the dimension of energy, A, through 

= ^7^- (3-13) 

We expect the gradient expansion to converge if the characteristic magnitude of the gauge 
fields is smaller than this scale. The temperature dependence of the parameter A is dis- 
played in figure 3. We observe a rather steep increase of A as a function of temperature: 
while at zero temperature, it acquires the relatively low value 53 MeV, at Tcfj = 1 GeV 
it has already reached the value 2.6 GeV. This is due to the fact that i^^jp is suppressed 
even more strongly than K^p^ at high temperature. Thus, we can conclude that nonzero 
temperature tends to improve the convergence of the gradient expansion. 
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Figure 3. The behavior of the scale A, defined in eq. (3.13), as a function of the effective temper- 
ature Toff. 

The definition of the scale A needs to be complemented by an evaluation of the magni- 
tude of the gauge fields in the order-six and order-eight operators. Ideally, we would want 
to compute numerically the relative size of these operators during a slow spinodal transi- 
tion. Such a calculation is beyond the scope of the present work, and to our knowledge is 
not available in the literature. The next-to-best thing is to try to estimate the approximate 
average magnitude of the fields during the transition. Below, we make three attempts at 
such an estimate based on different lines of reasoning. 

First, in ref. [27], the baryon asymmetry generated by a CP-violating operator of the 
type (4.1) was computed as a function of the speed of transition 1/tq, with tq being the 
quench time, that is, the time during which the mass parameter of the Higgs field flips 
its sign. The generated asymmetry turns out to be largest for itihtq < 20, where mn is 
the Higgs mass, corresponding to a characteristic momentum scale of /cq — 6-7 GeV. Our 
first rough estimate of the gauge field amplitude during the spinodal transition is that of 
a similar order of magnitude, say I-IO GeV. 

In the same paper [27], the SU(2) magnetic field was also computed as a function of 
time, as can be seen in figure 2 therein. Depending on the quench time and the time during 
the evolution at which we expect CP violation to be effective and the asymmetry to be 
generated, the magnitude of the field was found to be smaller than, or of the order of, the 
Higgs mass mn — 100 GeV, and decreasing with the quench time. Even for the slowest 
quench times shown there, the magnitude of the field is still considerably larger than the 
above estimated scale, A = 2.6 GeV, at Tcs = 1 GeV. 

Thirdly, we can use the results of ref. [28] to obtain a more quantitative estimate of 
the magnitude of the gauge fields. In that work, numerical simulations were performed 
within the somewhat simpler SU(2)-Higgs model, and the particle numbers of the SU(2) 
gauge field A^^ were computed. For the present discussion, we will assume that this gauge 
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field is representative of the field used in this paper. Accidentally, the computation 
was done in the unitary gauge suitable for comparison with our results here, but only for 
the fastest quench, tq = 0. The particle number can be converted to a correlator 

(4"(t.x)4(«,.)> = (l+ W_) J (3.14) 

where oj^ is the quasiparticle dispersion relation and the corresponding occupation 
number. (The superscript "T" here refers to the transverse modes of the gauge fields.) 
Taking for simplicity = ^Jk?' + and representing the occupation number by nj^ = 

lOe"^'^^!'*!/™'^ , as appropriate for the earliest time shown in the left panel of figure 7 in 
ref. [28], we get 

{A1{t,x)A]{t,x)) ~ 5''H,j{m GeYf. (3.15) 

Collecting the above three estimates for the magnitude of the gauge fields, and using 
the value of the scale A at Tgg = 1 GeV, we find that the ratio of the sizes of the typical 
type 6+2 operator and the type 4+2 operator is of the order 

/ (1-10^\ ' ^ ( < 100 GeV \ ' ^ ( 30 GeV \ ' ^ 

V 2.6 GeV ) -"-1^1^' ^ 2.6 GeV ) [2.6 GeV ) " 

(3.16) 

The third estimate is presumably the most reliable one, but it applies only to an instan- 
taneous quench. We conclude that although in the context of a spinodal transition, the 
gradient expansion may in some cases be poorly convergent, it would be desirable to do 
similar simulations for very slow quenches. 

3.3 Lepton sector 

As long as neutrinos can be treated as Dirac fermions, i.e., they possess no Majorana-type 
mass terms, a computation analogous to the above can be performed also in the lepton 
sector. Moreover, as there is no direct coupling between leptons and the strongly interacting 
sector, the calculation should be reliable all the way to zero temperature. Since the leading- 
order result is available from eq. (3.1), we merely need to plug in the appropriately modified 
values of the input parameters to evaluate the contribution of the lepton sector to CP 
violation. 

In the following, we take me = 511 keV, = 106 MeV and rrij- = 1.78 GeV for the 
charged lepton masses [9]. The absolute values of the neutrino masses are unknown, but 
using the fact that two differences of neutrino masses squared have been determined in 
neutrino oscillation experiments, Am|]^ = 7.50 x 10~^ eV^ and |Am|2l = 2.32 x 10~^ eV^, 
we can perform the calculation treating only the sum of neutrino masses as an unknown 
parameter. Finally, for the mixing angles in the lepton sector we take 9i2 = 0.59, 623 = 0.71 
and = 0.16 [9]. The CP-violating complex phase in the Pontecorvo-Maki-Nakagawa- 
Sakata (PMNS) mixing matrix is completely unknown; our approach is to fix it in order 
to maximize the lepton Jarlskog invariant, 5 = Tr/2. 

In figure 4 we show the zero-temperature coefficient Jk'^^ for the lepton sector, nor- 
malized to the same quantity in the quark sector, for a range of the neutrino mass sum 
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Figure 4. The ratio of the zero-temperature coefficient Jk"^^ in the Icpton sector relative to 
the quarks, for the whole allowed range of neutrino mass sum. The PMNS CP-violating phase is 
assumed maximal. 



up to 2 eV, roughly corresponding to the upper limit given by cosmological observations. 
The CP-violating coupling exhibits a clear power-law scaling with the neutrino mass sum, 
which is easy to understand. As long as the neutrino mass sum is much larger than the 
mass differences, neutrinos are almost degenerate with a common mass, mi,. Very much 
like in the quark sector, the momentum integral in eq. (3.2) can be evaluated approximately 
by factorizing out the masses of the heaviest contributing fermions, this time the charged 
leptons, leading to the simple expression 

Am^i(Am|2)^ 
CP ^ 807r2GFm»2- ^''■^'> 

The temperature dependence of the coefficient can be computed exactly the same way 
as for quarks. We find that the temperature suppression is even stronger for leptons and 
happens at much lower temperatures due to the fact that the lepton masses are so small. 
From the extremely small values of the lepton Jkqp^ ' conclude that their contribution 
to effective CP violation is negligible compared to the quark sector. For technical reasons 
related to the setup of our computational approach, this conclusion does not immediately 
apply if neutrinos have additional Majorana masses. However, we consider it likely that 
only extremely large Majorana masses would be able to compensate for the temperature 
suppression witnessed above. 

3.4 Checking the results 

Although a thorough cross-checking of the results prior to publication is a standard part of 
scientific work, we find it worthwhile to discuss it in some detail here. This is due to two 
reasons. First, we used a computer code that performs most of the calculations without 
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human intervention, and while this ehminates the possibihty of a trivial manipulation error, 
it may leave a feeling of a "black box" , the output of which must be analyzed particularly 
carefully. Second and more importantly, the description of our checking routines will 
provide the reader with additional insight into the structure of the result and its underlying 
symmetries. 

The leading order CP-violating effective action at zero temperature was originally 
computed in ref. [2] , where the operators not containing derivatives of the Higgs field were 
derived manually. Apart from verifying the consistency with this calculation, we performed 
two nontrivial cross-checks of our results. The first is based on itd-parity [2], an artificial 
discrete symmetry, under which the labels u, d of the fields as well as quark mass matrices 
are exchanged together with the replacement — )• W^f and — )■ —Z^. Such an 
operation can be compensated by a similarity transformation on the Dirac operator (2.4), 
so that it must leave the effective action invariant. On the formal level of eq. (2.7), it 
corresponds to a cyclic permutation of the terms in the trace. Therefore, verifying the ud- 
parity invariance of the result is a nontrivial check that the method of (covariant) symbols 
has been correctly implemented.^ 

Also the second, more nontrivial check can be formally understood as testing the 
implementation of the method of (covariant) symbols against a trace property that must 
be automatically satisfied. In contrast to cyclicity, this is the invariance of the trace under 
similarity transformations. As is obvious from eq. (2.3), both parity components of the 
effective action must be invariant under the transformation K — )• UKU~^ as long as U 
commutes with 75. A particularly interesting choice is U = (p/v, which commutes with 
everything inside K except for the covariant derivatives, which transform as UD^''^U~^ = 
D^''^ + (l)[D^''^ , = D^''^ — (p^. This shift of the covariant derivatives dramatically 
reduces the dependence of the intermediate result on the Higgs field, as it disappears 
from the offdiagonal components of the i^-operator (2.5) and thus no longer explicitly 
appears in eq. (2.7). While this similarity transformation only affects the dependence 
of the intermediate expressions on the Higgs field, it does constitute of a nontrivial test 
of the most complicated part of the calculation. As already remarked above, once the 
Higgs derivatives ip^ are discarded, much of the computation can be carried out manually. 
Why this is so can be seen from the sheer number of (/5-dependent terms appearing in the 
expansion coefficients (B.15). 

4 Applications 

4.1 Implications for baryogenesis 

The ultimate goal of our calculation was to pin down to what extent standard model CP 
violation can contribute to baryogenesis. Our result is particularly useful in this regard, 
since it not only provides a dimensionless measure of the strength of CP breaking (the 

^Note that while cychcity is a property of the full functional trace, it can be lost in the intermediate 
stages of the computation corresponding to eqs. (B.2) or (B.5), as the trace is performed only over the 
subspace of internal degrees of freedom in these expressions. 
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coefficients Cj), but also an effective, purely bosonic, action in terms of these coefficients 
and the corresponding operators. 

As we have seen, at temperatures above roughly 10 GeV, effective CP violation is 
suppressed by 10~^^ or more compared to the zero-temperature value, which leads us 
to conclude that in a high-temperature electroweak first-order phase transition, standard 
model CP violation is of negligible magnitude. Although the details of the calculation may 
change in the specific context of an expanding bubble wall interacting with a hot plasma, 
these conditions are unlikely to alter our conclusion, in particular at T ~ 100 GeV. At low 
temperatures, T ~ 1 GeV, the picture is less clear, as we will describe below. 

A temperature as low as 1 GeV for electroweak baryogenesis may be achieved if the uni- 
verse supercools way below the electroweak scale prior to electroweak symmetry breaking. 
This scenario of cold electroweak baryogenesis has been advocated for some time [16-19], 
and can be realized if the Higgs field is coupled to another scalar, which may [18, 29] or 
may not [30] be the inflaton, or if the first order electroweak transition takes place at a low 
temperature because the Higgs potential has a non-standard form [31, 32]. 

4.2 Matching to past simulations 

As mentioned above, the cold electroweak baryogenesis scenario has been studied in fully 
nonperturbative lattice simulations of bosonized effective theories, including SU(2) gauge 
fields, a Higgs scalar field and effective CP- and P-violating terms. The simplest such term, 

<5rcP-odd^ I -^^V'^VtrF'^'^F^., (4.1) 

where is the field strength tensor of the SU(2) group and (j) stands for the full Higgs 
doublet here, was considered in refs. [19, 27, 33], where it was found that the observed 
baryon asymmetry is reproduced for Scp — J [19]. Such an operator, however, does 
not appear in the standard model [14], at least in a gradient expansion of the fermion 
determinant. 

In ref. [34] simulations were performed with an operator of type 4+2, albeit of the type 
Oi , which — as we have demonstrated — is also absent from the standard model. Being 
based on a computation at zero temperature [1], the simulations did not take into account 
the i?i)-dependence of the (analogue of the) Cj's and hence diverged at (p = 0. To bypass this 
issue, a cutoff procedure was introduced [35], resulting in the asymmetry 

™B ,^_6 f A-50GeV\ ,^ 

-:.2.4xl0^x^exp^-^^^^j, (4.2) 

which upon correction for different normalizations of K^p^ here and in refs. [34, 35] roughly 
translates to (assuming the appropriate cutoff is around A = 50 GeV, as argued in ref. [35]) 

— ~ 1.5 X 10"^ X a, (4.3) 

where represents a generic coefficient of a 4+2 term. As the simulations using the full 
action (3.1) derived here are yet to be done, we make the very naive assumption that the 
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asymmetry generated from the Oj~-type term is comparable to the asymmetry produced 
by the P-even terms of eq. (3.1). The observed baryon asymmetry, eq. (1.3), can then be 
reproduced provided that Cj ~ 4 x 10~^, which roughly corresponds to T^g ~ 1 GeV. This 
means that successful baryogenesis is possible if the effective temperature at the time when 
baryon number violation is active is around or below 1 GeV, with an added contribution 
from each operator in the action (3.1). 

In a cold spinodal electroweak transition, the fermion fields are far from equilibrium, 
and we should thus strictly speaking redo the computation of the q's in the correct out- 
of-equilibrium state. As it is not clear how this can be achieved in practice, one might 
instead be interested in estimating an effective temperature, to gauge the strength of CP 
violation during the transition. This effective temperature corresponding to a given out- 
of-equilibrium state would then be whatever gives the same value for the q coefficients.^ 

In line with the discussion of convergence of the gradient expansion in section 3.2, one 
possible estimate of the momentum scale involved, and thus of the effective temperature, 
is obtained by translating the quench time tq into a frequency, /cq — Tq^- As mentioned 
there, we know from ref. [27] that in order to get a sizable asymmetry, the quench should 
not be slower than tquih — 20, thus leading to ~ 6-7 GeV. As we have seen that 
the effective CP violation is extremely sensitive to temperature, a more precise estimate is 
clearly necessary. Ultimately, the use of the bosonic effective action based on the gradient 
expansion should be replaced by first-principle simulations of the electroweak transition 
with dynamical fermions. In the following, we will discuss the future prospects for this. 

4.3 Future bosonic simulations 

Baryon number is C- and CP-odd, and P-even. Chern-Simons number is P- and CP-odd, 
but C-even. The chiral anomaly, which is a direct result of C- and P-breaking in the 
fermion-gauge interactions, relates the two through 

B{t) - B{0) = mcsit) - iVcs(O)], (4.4) 

and in order to generate a baryon asymmetry, they must therefore both become non-zero. 
This is only possible if all the discrete symmetries C, P and CP are broken. 

In a purely bosonic simulation, baryon number is not explicit, but assumed to follow 
the evolution of the Chern-Simons number. As a consequence, it is necessary to include 
the effects of both P and CP violation to generate non-zero asymmetry. Had the outcome 
of our computation of the effective action been CP violation in the P-breaking sector, the 
computed operators would be sufficient, as demonstrated in [34, 35]. But it has now been 
convincingly demonstrated [2, 20] that only CP violation in the C-breaking sector appears 
at leading order, and therefore bosonic simulations including only this would still give zero 
net Chern-Simons number. 

^Note that for very out-of-equilibrium states, this may be hard to achieve as the different a include hnear 
combinations of different moments of the momentum distribution. Reproducing all of the c; simultaneously 
by a thermal distribution is nontrivial. But considering a single Ci and a single operator, this can be 
envisaged. 
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Figure 5. The coefficient Kqp cio at Toff = 1 GcV for niud = as a function of nisc/mtb for 
mtb/Tctf = 1,2,4,8. Tlic successive values correspond to the plotted lines from the bottom to the 
top. 

The way out is to combine the P-even and C-odd operators from eq. (3.1) with the CP- 
conserving, C-odd and P-odd terms which originate from the left-handed gauge couplings 
to fermions. These should appear already at the fourth order of the gradient expansion 
and are naturally not proportional to the Jarlskog invariant. We have not computed these 
terms here, but the method of symbols is ideally suited for doing so, although complications 
may arise at this order through the Wess-Zumino-Witten term. 

As an aside, a similar complication arises in the two-Higgs doublet model. There, CP 
violation is provided at tree level through the complex couplings in the two-scalar potential, 
which break CP by breaking charge conjugation and preserving parity. Since CP violation 
is now located in the scalar potential rather than the fermion Yukawa couplings, one might 
expect to be able to generate an asymmetry using the bosonic fields only. However, parity 
is conserved, so the Chern-Simons number remains zero. One thus again has to include the 
C-odd and P-odd operators that allow for the generation of the Chern-Simons number. 

Simulations of this setup were recently carried out in ref. [36] , showing that an asym- 
metry is indeed created, but only when CP violation is combined with the C/P-violating 
operator (s). However, rather than explicitly computing the C- and P-breaking terms com- 
ing from the gauge couplings to fermions, the authors simply assume an effective operator 
of the type 

(^rc-odd/p-odd = i/^c/P {4>i<p2 - <^2(/>i) tr Ff.^F'"', (4.5) 
related to the chiral anomaly [37]. Here <j)i^2 are the two Higgs doublets. 

4.4 Future simulations including fermions: maximizing CP violation 

Ultimately, one would like to perform simulations including dynamical fermions in real time, 
along the way described in refs. [6, 7]. This bypasses the use of the gradient expansion, and 
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the fermions will automatically be in the correct out-of-equilibrium state. On the other 
hand, the numerical effort is huge, since the method requires stochastic averaging over 
fermion field ensembles for each member of a bosonic field ensemble. 

The observed baryon asymmetry is n^/n^ ~ 10~^, and such a signal is very hard to see 
in numerical Monte-Carlo simulations. In practice, what is done is to simulate the dynamics 
with a value for (5cp (of whatever bosonic operator) much larger than the physical one, find 
the regime where the baryon asymmetry is linearly dependent on 5cP) and interpolate the 
result to the physical value. For purely bosonic simulations, one can in principle choose 5cp 
to have any value. When including fermions, the CP violation is however encoded in the 
actual CKM matrix, the structure of which does not allow arbitrarily large CP violation. 
What we can do is try to maximize the baryon asymmetry by modifying mixing angles, 
the CP-violating phase and the fermion masses. Although we do not a priori know how 
the asymmetry depends on these parameters, we can get a good estimate of the optimal 
point in parameter space by maximizing Kqp^ 

First of all, CP violation is generically proportional to the Jarlskog invariant, which 
takes its maximum for the (unphysical, but in principle allowed) mixing angles sin0i2 = 
l/\/2, sin 013 = l/\/3, sin023 = l/\/2 and the complex CP-violating phase 5 = tt/2. The 
resulting value J = l/(6\/3) ~ 0.096 is larger than the physical value by the factor of 
3.3 X 10^. 

Second, the CP-violating effective couplings depend in a particular, and rather non- 
trivial, way on the quark masses. The overall factor A is simply maximized for large mass 
separations, but the complete expression is more complicated. We start by simplifying the 
problem by taking the n-type and d-type quark masses to be equal, m„ = = m^^, 
rric = rUg = rUsc and mt = nib = mtb- In addition, we will take the limit mud since 
the actual value of this mass plays very little role; gular and does not change sig- 

nificantly in this limit. Consequently, the coefficients k'^Ci depend on three independent 
scales: m^c, ^Mb and the effective temperature T^s. 

By simple dimensional analysis, be written as l/^Jj times a function of two 

dimensionless ratios, which can be chosen conveniently as irisc/mtb and mtb/T^s- While 
the former characterizes the fermion mass hierarchy, the latter compares the overall mass 
scale to the temperature. Figure 5 displays the largest of the coefficients, 
function of the mass hierarchy nisc/mtb for several values of mtb/T^s- While the plot applies 
directly to the case of Tgfj = 1 GeV, the values of the coupling at other temperatures can 
be obtained by mere rescaling. 

Varying the three available parameters, Tefj, rUsc and mtb, creates a lot of space for 
enhancement of the signal for CP violation. For instance, by choosing T^s = 1 GeV and 
mtb = 4 GeV, figure 5 tells us that we should tune rUsc so that msc/mtb ~ 0.6. The 
coefficient K^p^cio then reaches the value of roughly 10~^, which is to be compared to 
its value at the physical point, 3 x 10~^. While this is merely a concrete example of the 
dependence of the effective coupling on the parameters as displayed in figure 5, additional 
constraints due to the setup of the numerical simulations need to be taken into account. 
The fermion lattice simulations in refs. [6, 7] acquire much larger statistical uncertainties 
for large fermion masses, so we would like to keep nisc < mritb ^ 10 GeV. On the other 
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hand, we also need to resolve the physical gauge and Higgs boson masses on the lattice, so 
the fermions should not be too light, either. 

By combining the above concrete choice of fermion masses with the optimization of 
the Jarlskog invariant, we are able to increase the expected asymmetry by a factor of at 
least 10^. Conversely, if one performed simulations including fermions with these optimized 
parameters, then finding a net asymmetry of the order of 10~^ should correspond to the ob- 
served asymmetry at the physical point. Note that this is precisely the order of magnitude 
of the asymmetry seen in bosonic simulations (see for instance ref. [33]). In other words, 
although the scope to tune CP violation is limited to only shifting masses and mixings 
rather than freely tuning 6cp, the numerical signal can still be made large enough to be 
seen in a real-time Monte Carlo lattice simulation with numerically manageable statistics. 

5 Conclusions 

Integrating out fermions from the standard model gives rise to effective CP-violating op- 
erators for the gauge and Higgs fields. In the present paper, we have shown how one can 
obtain such operators systematically within a covariant gradient expansion, at both zero 
and nonzero temperature. We have computed the temperature dependence of the eff'ec- 
tive couplings at the leading nontrivial, sixth order in the gradient expansion, and found 
that to this order, there are no operators that simultaneously violate both the CP and P 
symmetries. Although detailed numerical simulations using the resulting bosonic effective 
action are yet to be performed, a simple estimate allowed us to conclude that standard 
model CP violation may be the source of the cosmological baryon asymmetry provided 
that the electroweak phase transition happens at low enough temperature, of the order of 
1 GeV. This implies that while "hot" electroweak baryogenesis with the standard model 
as the only source of CP violation is firmly ruled out (as has been known for many years), 
other cosmological scenarios such as cold electroweak baryogenesis might provide a viable 
solution to the baryon asymmetry problem. 

With regards to numerical simulations of baryogenesis, we made several observations 
that will partially direct our future efforts. First, the contribution of the lepton sector to CP 
violation can be safely neglected as long as neutrinos are of the Dirac type. The discussion 
of Majorana neutrinos goes beyond the scope of the present paper and will require a 
proper extension of the formalism used here, but we suspect that our qualitative findings 
will not be altered. Second, numerical bosonic simulations determine the generation of 
the gauge field Chern-Simons number, which is related to baryon number via the chiral 
anomaly. As the Chern-Simons number is odd under parity, and all of the operators we 
have found in the leading order CP-violating effective action preserve parity, the knowledge 
of simultaneously P-odd and C-odd operators is required. These arise from the purely left- 
handed gauge interactions of fermions and are expected to appear already at the fourth 
order of the gradient expansion. 

As a final technical point, we have estimated the optimal values of fermion masses 
and mixing angles to give the maximal baryon asymmetry. This is important for future 
simulations of baryogenesis with dynamical fermions and the actual CKM matrix structure, 
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since we then do not have the freedom to make the asymmetry arbitrarily large (and 
numerically detectable) by choosing a large 6cp freely. 
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A Notation and conventions 



Throughout the paper, we use the natural units, in which the Planck and Boltzmann con- 
stants as well as the speed of light are equal to one. We also work in Euclidean space; 
Euclidean spacetime vectors are denoted by lowercase italic letters and their spatial parts 
in boldface, for instance v = {vo,v). Spacetime indices are represented by lowercase Greek 
letters; repeated indices are summed over unless explicitly stated otherwise. In the Mat- 
subara formalism, the temporal component of momentum is discrete. For fermions, it takes 
the values Un ^ (2n + 1)ttT, where T is the temperature. 

In order to simplify the notation, we adopt the following conventions for integration 
symbols. The integration over spacetime coordinate x is denoted as = f d'^x at zero tem- 
perature, where d is the dimension of the spacetime. At nonzero temperature, it has to be 
modified to = f^^'^ dxo J d'^~^x. Momentum integration is defined as = / d'^p/{27r)'^ 
at zero temperature, and f = T^^^^f d'^~^p/{2Tr)'^~^ at nonzero temperature. Here, the 
sum is over the Matsubara frequencies, pQ = Un, as defined above. 

A.l Standard model 

In order to simplify the lengthy expressions we have to deal with, it is crucial to introduce a 
suitable notation. We follow ref. [2] and absorb the gauge couplings into the redefinition of 
the gauge fields. Specifically, we trade the physical fields Wj^ (charged weak intermediate 
bosons), (neutral weak intermediate boson), (hypercharge gauge field), and Gafi 
(gluon color-octet field) for 

= -^W^, = ^——Za, = g'B,„ = —XaGan- (A.l) 

^ ^ ^' ^ 2cOs6'w M y ^ a a^l \ ) 

Here, as usual, g,g',gs denote, respectively, the weak isospin, hypercharge, and strong 
coupling constants, while is the Weinberg angle and Aq are the Gell-Mann matrices. 
With this notation, the chiral components of the (Euclidean) quark Dirac operator (2.2) 
in the unitary gauge become 



W- Dd,,,-Z^ + Gj' \ fMrf 



^^^^^ z.,, + gJ' --Mo 



(A.2) 
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where the matrix structure corresponds to isospin space. The hypercharge covariant deriva- 
tives used here are defined as -Du,/^ = + {2/3)B^ and -D^,^ = — {1/3)B^. In the 
two-dimensional space of eletrically neutral electroweak gauge fields, it is common to work 
with the basis of the field and the photon field, or with that of the neutral isospin gauge 
field and the hypercharge field B^. We choose here a "mixed" basis, composed of and 
B^ . This is the reason why Z^ does not appear in the right-handed component of the Dirac 
operator, and why the couplings of the hypercharge field in the covariant derivatives 
equal the electric charge of the fermions rather than their actual hypercharge. 

For the sake of completeness, we also write down the analogous expressions in the 
lepton sector, 

""'^^^y w- dJ-z,)^ ---^^ imJ ^^^^ 

where we denoted -De,/^ = dfi — B^ and assumed that the neutrinos have only a Dirac-type 
mass matrix, My. 



B Method of symbols 

The method of symbols is a systematic device to calculate traces of differential operators. 
While its basic form has been known for a long time [38] and its covariant version was 
introduced more than a decade ago [39], the generalization to finite temperature has only 
been worked out recently [25]. Here we provide an overview of the method, following closely 
the notation of ref. [2]. 

Our task is to calculate the trace of a differential operator f{D,M) where is a 
covariant derivative and M{x) a (local) background field. Both of them can in principle 
have some additional matrix structure (internal degrees of freedom). In order to calculate 
the trace we need to choose a basis. It is convenient to choose a direct product basis, 
composed of independent bases in the space of test functions of the coordinate x, denoted 
as Xn{x), and in the internal space. Provided the basis Xnix) is orthonormal, the full trace 
becomes 

Tvf{D,M) = Y, / tv[xUx)f{D,M)xn{x)], (B.l) 

where "tr" refers to a trace over the internal space only. Choosing in particular the basis 
of plane waves, Xp{x) = e'^'^, we replace the sum over n with an integral over p and use 
the operator identity e~'P'^De'P'^ = D + \p to obtain 

Tr f{D,M) = [ tv[f{D + ip,M)l], (B.2) 

J x,p 

where 1 is a unit matrix in the internal space. This formula constitutes the main result of 
the method of symbols. 
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While the above expression is fully general, its obvious disadvantage is the lack of 
manifest gauge covariance due to the appearance of the "free" covariant derivative oper- 
ators [2]. In order to remedy this, let us first introduce an arbitrary vector operator 
acting in x-space and rewrite eq. (B.2) as 

TTf{D,M)= f tr[e-'^-V--^e-^P-^/p,M)e'PV^-^IL], (B.3) 

J x,p 

where we used the shorthand notation = d/dp^. Note that adding the exponential on 
the far right does not change anything since it acts on the unit matrix, while the two extra 
exponentials before f{D,M) combine to a unity. Therefore, this is an identity that holds 
regardless of the temperature. 

Next, we define a barred symbol for an operator v through 

V = e'^-^e-'P-^ V e^P-^e"'--^. (B.4) 

As the last step, we expand the first exponential e~'~'^ in powers of the derivative and note 
that its spatial part produces zero, as it gives an integral of a total derivative in momentum 
space. (This observation relies on the fact that acts in coordinate space only and thus 
commutes with V^.) The final result for the trace of the differential operator in the method 
of covariant symbols therefore reads 

Trf{D,M) = J2^-^ tr[(HoVor/(AAf)l]. (B.5) 

n=0 ™- -^^'P 

Note that at zero temperature only the n = term of the series survives, and in addition 
that we have complete freedom in the choice of the operator H^. It is, however, advan- 
tageous to simply set = D^, since the "free" derivative in + ip^ is then cancelled 
and Dfj^ becomes a multiplicative operator, containing covariant derivatives only inside 
commutators. Concretely, the expansions of the covariant derivative and the scalar field 
read 

oo 

D, = e'^-^(Z)^ + ip^)e-'^-^ = il^M + E (n + l)(n - 1)! ^"""""^^"^ " " " 



^— ' T7 .I 



(B.6) 



n=l 



= M + i[D„, M]V„ - ^[D^, [Dp, M]]V^Vp + • • • , 

where M«,...„„ = [D^,, [■ • • M] • • • ]] and F«,...„„^ = [D^,, [■ ■ ■ [Da^,D^] •••]]. This 
concludes the derivation of eq. (5) presented in ref. [20], and agrees with eq. (2.40) derived 
independently in ref. [25]. 
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B.l Application to the standard model 

In order to be able to apply the method of covariant symbols to the standard model, we 
have to respect the isospin structure of its Dhac operator. We therefore choose 





Df,^ 



(B.7) 



in isospin space. Should the terms including the gluon field be desired as well, we would 
simply replace D^'"^ with D^''^ + in the above. As the gluon field commutes with all 
other fields but the covariant derivative itself, this shift would only affect the field strength 
tensors Fa'^..oi„^ defined below, giving them a separate contribution in color space. The 
barred electroweak fields, defined by eq. (B.4), now take the form 

n=l 

oo 

li. 

n=l 

oo 

Z^i = e'^ ' '^Z^e~'^ ' = + ^ —Zai-a„fj.^ai ■ ■ ■ Va„, (B.8) 

n=l 
oo .„ 

^^l = e (Pf,e = </2/i + 2^ — V^ai-anM^ai " " " 

n=l 

where F^'^..^^^, = [Daf, [■ ■ ■ [Dalt, D]^''^] •••]], Zq^.^^^^ = [da, ,[■■■ [da„,Z^] •••]], and anal- 
ogously for •^al■■■ar,^l■ Furthermore, we have defined here VF+02-a„M = ~ 
and Wa,a2-a,,t, = Di,W-^-a^^ - ^a2...a„M-^m- ^ote that the above com- 
pact expression for Wj^ follows from the generalization of the Baker-Campbell-Hausdorff 
formula, 

e^Ae-^ = exp([^, •])A, (B.9) 
where A] = AX - XB. Finally, we need the expression for (p, 



oo \ 
n=l ' ^ 



where the structure of the J7-tensor is determined by the partitions of the set of indices 

Q?! , . . . , Cln) 

partitions i 

Ui{ftj}j={"n}" 
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for example 



^al3 ='Pal3 + V^a^(3, 

^ajS-yS =y^afS'y5 + 9^/3 9^7 '/'^ (B.12) 

In order to be able to evaluate traces of the type (2.7), we still need to know the 
propagator in the method of covariant symbols, cf. eq. (2.6), 

k!d = ^Mu^mI, - {Jl)u4 ± t\^uA + (B.13) 

This is expanded up to the desired order in the covariant gradient expansion, -/V~^ = 
^ud ~ {Ui''^ + f^2''^ + ■■■)' ™d the series is inverted so that the propagator reads 

= N^4 + NuAUi''' + U^''')Nu4 + Nu4U'i'''n^4Ui'''n^4 + ^^^ terms. (B.14) 

The expansion coefhcients of the inverse propagator read, up to fourth order in the gradient 
expansion, 

- ± hz^p,h>ycNtif - -AhKp, + yc.p,)fi>yc.vp+ (b.i5) 

jju,d ^ _ Gu^ij^falSyS + I'falSyfS + ^H^apVyS + '^''faP'fy'fS + \'^aV P'^y'^d)^ (N y^ 5+ 

- (lF„";f ± (If;;'^ + 

where we have introduced the shorthand notation Gu,d = ^ud ~ P^' This ensures that all 
momentum integrals resulting from the method of covariant symbols are of the type (2.9). 
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C Momentum integrals in the 4+2 sector at finite temperature 

Upon diagonalization of the quark mass matrices, the momentum integrals needed in the 
evaluation of the 4+2 type operators in the CP-violating effective action take the form 



"/^ 3— /c . V — ^ / 
A, 1,2, 2 ~ i Ini y ^ Jijmn I 



i,j,m,n 













(p2 + 








,mY{P^+m\^Y 








-k 




(p2 + 








.mYiP^+^Xuf 






ipyipl)'- 


-k 




(p2 + 




^\j){p'^ + 


ml. 


.mfip'^ + r^lnf 



/''^-'^-ilm V J / ' ^^0^ (ru 

i,j,m,n P 

A,'l,3,2 = i-l^™ ^ Jijmn J 
i,j,m,n P 

see eq. (2.9). Here we have defined Jijmn ^ ^ij^jni^mnVni^ , for which eq. (1.2) gives 
Im Jijmn = J^im^jn- The notation used for the quark masses should be obvious; here and 
in the following, the indices i,j,m,n always assume the values 1,2,3. 

To simplify the above expressions, we first use an identity given in eq. (11-2) of ref. [2], 

\ " ^im^jn _ 

which allows us to rewrite i^f f 2 2 form 



ntA^uJiiJ^) f y^w" . (c,3) 



Similarly, we obtain for the sum /{^'f 23 + -^1 1^32 i'^^^Y ^^i^ combination appears in the 
effective action) 

1 d 



jkA-k , jkA-k _. J Sr^ [( 2^^k( 2^i-k 

-'1,1,2,3 -'1,1,3,2 —l'^ timtjn / IP j \Pq) 



i,j,m,n 



2 9p2 

^ ■ (C.4) 

X 



(p2 + m2^J(i52 + m|^.)(l'^ + ml,niV{p'^ + 'ml ^2 ' 



which can be cast in the form 



A,i,2,3 + A,i,3,2 - ^ 1^115^2 J y^n:=i(P^+-i)" 4 '^-^'^'^ • ^^-'^ 

From these expressions, we observe that we only need to evaluate the class of integrals 

IpWU^P^ + ml 



as well as their mass derivatives 

/ 6 

r 



\q-- 
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for k = -1,0,1,2,3. 

It would clearly be desirable to evaluate the above integrals in as symmetric a fashion as 
possible, and preferably in a way, where it is not necessary to specify the number of flavors 
A^f (A^f = 6 in our case). The calculation is then somewhat simplified by the introduction 
of a new set of integrals, 



'(0) _ 




k — 


I- 






'W - 


1 


k — 


iVf 


-m - 





2\k 



/^2\k 



1 ^ f (p2)fe 



TT\-'J = - \ ^ f 



(C.8) 



" iVf(iVf - 1) • • • (iVf - n + 1) ^^^^^^^ I n,^{,,...,,„}(P^ + m^g) ' 

where the angular brackets indicate that the sum is taken over sets of mutually different 
indices only. The absence of the frequency variable in the numerator makes these integrals 
particularly easy to compute. To evaluate them, we next define 



1 . Er, 



Eg = Jp^ +ml, ag = — tanh— ^, Pq = Yi 2T' (^-9) 

V 2Eq 2T^s Ug'^gim^g - ^g') 

as well as the new mass parameters 



1 {q,q) 

' " ~ Ni{N{ - 1) • • • (TVf - n + 1) ^ ■ 

{qi,---,qn) 



(C.IO) 



Using these quantities, we straightforwardly arrive at the expressions 
which are painless to integrate numerically. Fully analytic expressions can in addition be 



(C.U) 
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obtained at zero temperature, 



- 4^5/2 rOMh^ z^'^i A 



\k 



(1) _ (-1) r(3/2 + fc) ^ 2(l+fc)(^^2 _ ^2^o 1 

\ ' q 

(2) _ (-1)1+^^ r(3/2 + fe) ^ ^2(l+fc) /„ 4 9^2..2 , ,.4xo , 
-^^57^ r(2 + fc) 2^""? '("^g - ^"IgAii + j/:;^ log — 



(C.12) 



^ ^ q 

where the (arbitrary) parameter A is introduced to make the argument of the logarithm 
dimensionless but drops from the result upon summation over q. As evident from these 
results, all integrals defined above are both ultraviolet and infrared convergent in four 
dimensions. 

As the last step of our derivation, we then note that the integrals Ik can be expressed 

(n) 

in terms of the Hj^ via 

/_i =24Mf//i°; + QMlH^^l - MlH^^l + AMIh^^^ + H^^ + 24Mf (i^^^^ - i/^) 

- 4M|ifi\^ - AH^^^ - 12Mf{3MiH^^} - ^ + 2H^'^ - H^^l) 

- QMl{Hf^ - 2H^^^ + ifi'^) + QHf^ + AMI [2MIh^^} + iff ^ 

- 6M|(ii(°) - - 3i/« + S/if - H^^l] - 4//f + Fit (C.13) 
Jo = - 6Mfif(°) - M|iii°) + 3M|iif ) - /if - 6Mf{H[''^ - h!,'^) - SM^hI,'^ 





+ 3i/« - 


3Aff (-2M|i7j°^ 


h 


=2MfHf^ - 




h 


= - MlHf^ 




h 







A3) 



+ 2Mf{Hi'^-Hi'^)-2H^'^+Hi'\ 



From here, the desired integrals iff 2 2' -^1123 ^^'^ -^1132 reconstructed using eqs. (C.3) 
and (C.5). Finally, we will for completeness reproduce the explicit expressions of the 
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coefficients q, appearing in eqs. (3.3)-(3.5), in terms of tlie 1^ integrals, 

/3 _ 35/1 + 70/2 + 27/3 _ 7/2 + 3/3 

(•••)|t=o' (•••)It=o ' (•••)It=o' 

35/1 + 42/2 + 15/3 35/1 + 42/2 + 23/3 35/1 + 70/9 + 19I3 

C4 = 7 71 , C5 = -. , C6 



(•••)|t=o ' ^ (•••)It=o ' ° (• 



|T=0 



_ 35Ji + 28/2 + 9/3 _ 3/0 + 7/1 + 5/2 + /3 _ 35/1 + 14/2 + 11/3 

(•••)|t=o ' (•••)It=o ' (•••)lr=o 

105/0 + 105/1 + 119/2 - 9/3 105/1 + 42/2 + /3 ,^ , 

Clo = -f , cii = , (C.14) 

(•••j|T=0 (•••j|T=0 

8(7/2 - /3) 4(35/1 + 14/2 - 5/3) 

C12 = —7^ rr. ; C13 — 



35(/o + 3/1 + 3/2 + /3)|t=o 35(/o + 3/1 + 3/2 + /3)|t=o 

Here, the notation in the first 11 coefficients indicates that they are normalized by the 
respective zero-temperature results. The normalization of Ci2,ci3 has been chosen in such 
a way that the prefactors of the corresponding terms in eq. (3.4) are equal to one. 

D Momentum integral in the 6+2 sector 

Similarly to eq. (C.l) of the previous appendix, the momentum integral needed for the 
evaluation of the CP-violating operators of the 6+2 type at zero temperature reads 



^1,1,1,1,2,2 = ilm ^ Kijklmn / {p 



1 



l'm,n (D.l) 



(p2 +m2 „)2(p2 _^m2 „)2' 

where Kijkimn = ^ij^jk^^kiVi^VmnVj^^ , analogously to the Jarlskog invariant defined in 
eq. (1.2). In order to evaluate this sum, we first try to simplify the expression for the 
/^-tensor. 

Suppose first that i = k. Thanks to the unitarity of the CKM matrix, one finds 
ImKijiimn = \Vij\'^ImJii^n = J^imGlnlVijl'^, thus proving that for i = = m, Kijkimn is 
real. We conclude from here that its imaginary part can be written as a sum of terms 
corresponding to three possible cases, in which two of the indices i,k,m are equal, and one 
case, in which they are all different, 

Irn/^ijfcimn — J^ki^km^ln\^ij\ ~^ J^mk^mi^nj\^kl\ ~^ J^im^ik^jl\^mn\ 

+ \e,km\Mv,,v-,'VkiVi;;!VmnV-'). 

In the last term, we can use the unitarity of the CKM matrix to eliminate the index m, 
Vi~^Vmn = Sin " Vii^Vin " ^fc^V'fcn, leading us to 



llll Kijkijjin — ^ki^km^lnl^ijl ~^ ^mk^mi^-njl^kll ~\~ Sim^ik^jl\^mn\ 



1 

J LJh,lllLlb ret - l\,lll,-l,IL\ ' tj \ • " IILrC-lltflLJ \ ' h^t \ • " L 1 1 L - Ihl, — J I \ ' IILIL \ (D 3^ 

H" \^ikm\^ik^jl{_^ln |^n| ) | ^^fcm | ^ifc^jT^ | V/^^/ 1 . 
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This form can be further simphfied using the trivial identity {cikml^ik = ^ikm the fact 
that Kijkimn is invariant under cyclic permutations of the pairs of indices ij, kl and mn, 
allowing us to write the tensor if in a maximally symmetric form 

— \in. KijjflYnn — ^ki^km(-ln\Vij\ ~l~ ^»nfc j ^nj | | ~^ ^im(-ik^jl\^mn\ 



J 



+ ^eikm [ein {Snj + \Vij\^ - \ Vkj\^) + enj {Sji + I Vfc; 1 2 - I T/^; 1 2) (D.4) 

Plugging this expression back to eq. (D.l) and performing some straightforward albeit 
tedious manipulations, we arrive at the final expression for the momentum integral in the 
6+2 sector, 

which is simple to evaluate numerically. 
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